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Abstract 

We combine neutron scattering (INS) data and NMR/NQR nuclear spin lattice relaxation rate 
(1/Ti) data to deduce the existence of a new contribution to the magnetic response ~x"{q,uj) in 
cuprate superconductors. This contribution, which has yet to be observed with INS, is shown to 
embody the magnetic pseudogap effects. As such, it explains the long-standing puzzle of pseudogap 
effects missing from cuprate INS data, dominated by stripe fluctuations, for x"(<f, to) at low energies. 
For _Lai.86<SVo.i4C'u04 and YBa2Cu^O^^, the new term is the chief contributor to 1/Ti for T ^> T c . 
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Cuprate superconductors are well-known for their unusual normal metallic state proper- 
ties. Prominent among these is an extended 'pseudogap' regime located above the dome of 
superconducting transition temperatures T c (see e.g. a typical phase diagram [1]). Theorists 
have attempted to connect features in the normal state phase diagram with the phenomenon 
of high temperature superconductivity itself. There are three main approaches. The first 
focuses on the disappearance of the Mott-Hubbard antiferromagnet, which is transformed 
into a valence bond state where mobile holes are naturally paired [2|. The second concen- 
trates on 'stripe' correlations as providing the environment needed for superconductivity j^]. 
The third hypothesizes that the cross-over into the pseudogap regime is actually a phase 
transition to a state with concealed long range order (4]. Experimentally, we have been 
able to use inelastic neutron scattering (INS) as a function of temperature, composition and 
magnetic field to map stripe order and fluctuations . In the joint NMR and INS analysis 
presented here we identify a novel low-frequency signal which is correlated with pseudogap 
formation. 



Beginning as an NMR effect [SHlO], the pseudogap has been observed as a genuine charge- 
energy gap llJ-JJal. with excitations known as 'arc fermions', characterized in detail through 

! ~™ nn 

recent ARPES studies [14|, [15[. Thermally induced changes of the Fermi surface, with 
concomitant behavior of arc fermion excitations, are clearly related to magnetic manifesta- 
tions of the pseudogap. However, INS data for x"(g, w) show only indirect manifestations 
of the pseudogap [la, Meanwhile, systems such as Lai. 8 6Sr .i4CuO4 (LSCO) [a] and 
YBa 2 Cu306.5 (YBC06.5) j3| yield data for x"(g, w) that consists at low frequencies of in- 
commensurate, antiferromagnetically correlated peaks whose intensity exhibits cu/T scaling 



181 ] from T~ 60K up to room temperature. Pseudogap effects are totally absent from such 



data. Moreover, nuclear spin-lattice relaxation rates (1/Ti) for these systems 19M21| are 
inconsistent with extrapolation of the INS results to NMR frequencies. 

In this Letter we present a new, joint analysis of INS and NMR (T\) data for the systems 
mentioned above, in which we deduce the existence of a pseudogap fluctuation term Xp(<f> u )i 
which has not been identified by INS up to now. Thus, we write x"(<f, co) = ^"((f, u;) a + 
Xp{q,u)), where x'liQi^a is the INS-measured term with incommensurate peaks along the 
a axis. The term x'p(Qi u ) introduced here, which is nonzero in the fluctuating stripe phase, 
will be modeled below to interpret the 7\ data. Thus, not only does the strongly evidenced 
occurrence of such a term clearly explain the hitherto baffling omission of a pseudogap effect 
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from data for u;)p, |7|], it also accounts for the disparate behavior of T\ for the planar 



63 Cu and 17 nuclear spins in these systems [20, |21|. We also show that the thermal and 
q-space behavior of Xp(<f, is such that it could have easily been missed up to now by INS 
experiments on these systems. In sum, the present analysis addresses a major deficiency in 
our understanding of the anomalous normal-state physics of cuprates and will be of wide 
interest to theorists and experimentalists alike. 



For the nuclear re 
by Uldry and Meier 



axation analysis we employ the formulation of the T\ process pioneered 



221 ]. in which the relaxation rates are written 



63 Ti 



A 2 ah + AB 2 + 2B 2 {AK 2 + K 3a + K 



36 J 



+ AA ab B(K la + K lb )}r eff 



17 T lc 



[C 2 a + C 2 }(2 + K la + K lb )r ( 



eff 



(1) 

(2) 



for bi Cu and O, the two nuclear species of interest. In these equations A ab , B, C a and 



C b are hyperfine tensor components in units of Gauss per unit of spin 
K n represent the normalized dynamical spin-spin correlation functions 



2 011 . T he quantities 
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24]. Thus, K n 



A(Si ■ Sj), where Si and Sj are n neighbor spins. The additional subscript a, b indicates, 



for n = 1 and 3, that the bond axis (i.e. 



fluctuation-dissipation formulation of T± 
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ies along the a or b crystalline axis. Using the 



261 ] . we express the K^s in terms of x"(<7> w )a, 



taking a as the discommensuration axis. The K n 's (n = 1,2,3) for the first three neighbor 
pairs in the CuO"2 plane may then be written 

In dq g n a,b(q) [x"(q, u)J u]^o 



K 



na,b 



(3) 



I n dq\x!'(,Q,u)aMu-*o 

where gi a>b = cos(q atb a); g 2 = cos(q a a)cos(q b a); and g 3a ,b = cos(2q atb a) in an obvious notation. 
Note that K 2 is independent of the discommensuration axis. It is clear that \K n \ < 1 . Eq.flT]) 
and (|2J) also employ the key parameter 

k B T 



Teff(T) 



\1 B Jn 



f dq 


x"(q,u) a 


In 


UJ 



(4) 



J u;->-0 



proportional to the 'local susceptibility', which acts as a correlation time that includes the 
particle statistics of the relevant carriers. r e ff may also be estimated directly from INS data. 
Thus, we shall proceed by comparing the latter values of r e ff with those obtained from Ti 
data, using reasonable estimates of the other parameters in Eqs.Q and (j2J). 
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We begin by applying these equations to LSCO, using INS data for x"(Qj ^Opi designated 
x'j(q,uj) a , in Eq. (j5J) to estimate r e //, which we denote r e j. The resulting values |27|, plotted 
as solid triangles in Fig. 1(a), exhibit approximate u/T scaling (i.e. T e ff = constant) above 
the stripe onset temperature Tg tr ~ 50K [18 1, in accord with the fluctuating stripe domains 



model of Zaanen et al. 
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29]. Below Tstri r e/(T) drops essentially linearly toward zero. 
We note that T c ~ 35K for this sample jfj|. 

The latter results are to be compared with estimates of r e //(T) derived from T\ data 
20| using Eq.(l) and (2). To do this, we employ values of the HF constants derived from 
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shift data [30j and values of the K na ^s calculated with Eq.flH]), using a two-peak form factor 
fitted to INS data [6]. The K nIa ^s so determined vary only gradually with temperature. 
As an example we mention values at T = 100K: Ku a = -0.75; K^ib = -0.97; K 2 i = 0.73; 
K^ia = 0.16; Kzit, = 0.90. One notes a sharp distinction between the a and 6-axis values 
as expected. For the 17 O, this could result in the prediction of two widely different values 
of l/ 17 Ti c for oxygen sites with Cu-0 bonds aligned with the a and b axes, respectively. 
Since only a single rate was observed {20]] , it is presumed either that a flip-flop mechanism is 
present to maintain a single 17 O nuclear spin temperature, or the stripe domain boundaries 
are fluctuating, so that each site automatically averages the two rates to yield a composite 
rate (given by Eq.@). Results so obtained are plotted in Fig. 1(a) as blue squares and red 
circles for the 63 Cu and 17 nuclear spins, respectively. For comparison, the value of r e jj 
deduced for YBC07 by Uldry and Meier 22| is shown as a dashed line, reflecting the fact 
that the T\ process in LSCO is substantially stronger than that for YBC07. 

Regarding the estimated values of r e // for LSCO in Fig. 1(a), note that there is approx- 
imate Korringa-like behavior for all three data sets below T = T$t r , where the agreement 
is good considering that there are no adjustable parameters. The slight disparity in mag- 
nitudes is attributed to HF constant errors or differences in the widths of incommensurate 
peaks (i.e., of the K n 's) among samples (see the YBC06.5 case below). Since the curves for 
T e ff rise high above those for r e j at T > Tstr, the T\ data give clear evidence for an additional 
term in x"(<f, a;) as stated above. We emphasize that there is no other realistic possibility 
among the well-documented T\ mechanisms in solids. Also noteworthy is the drastic differ- 
ence between the r e ff curves derived from the 63 Cu and 17 T\ data. Values of r e fj for these 
two measurements must actually be the same, so that the correlation properties of the new 
term x'piMi 1 ^) are evidently rather different from those of X/(^i a; )a- 
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Our next step is to take explicit account of the two terms in x"(q> u) and rewrite Eqs.(l) 
and (2) as 

1 _ 1 

637T 6371 v ' 

lc -«■ 1 



1/c 



and 
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2 
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A^ 6 + 4^(1 + 2K 2P + K 3P ) + 8A ab BK 1P 



, 2 



l7;/ , ( : T^ + ^[Ca + C b 2 ]((l + ^ipKp, (6) 

where !T n p and r e p are defined using x P (g, a;) in Eqs.(3) and (4), respectively. In Eqs.(5) 
and (6) the Ti/ C 's are calculated with Eqs.(l) and (2), respectively, using r e // = r e j, and 
= K n j. We may now use Eqs.(5) and (6) to extract estimates of r eP from data for 
both l/ 17 Ti c and l/ 63 T lc . However, it is necessary to model XpiQj^) i n order to make 
systematic estimates of the K nP . For this purpose we follow Aeppli, et al. Jfj], taking a 
squared Lorentzian form of unit amplitude Xp(Qj / Xp(.P ea k) = qt/ilw + 1x + <?y) 2 > from 
which the K nP s follow via Eq.([3]). Since 17 T" lc (Eq.(jS])) varies rapidly with K iP , while 63 T lc 
(Eq.(5)) is more weakly dependent on the K nP s, the width parameter q w may be varied 
with temperature to bring the r eP s into coincidence. Results of this procedure for LSCO 
are presented in Fig. 1(b), with the corresponding values of K nP shown in the inset. 

The squared Lorentzian form for XpiQ^) centered on (tt,tt) gives a satisfactory account 
of the data, where we have taken, e.g., K\ P (T) = -0.8lexp[— (T — 50)/600]. The width 
parameter varies between q w ~ 0.6 and ~ 1.3 (units of a -1 ) for 50K < T < 300i^. q w is 
therefore similar to the displacement of the incommensurate peaks (~ 0.77) in LSCO. 

Next, we consider the case of YBC06.5. While LSCO has only a weak pseudogap, 



YBC06.5 has stood from the earliest days as a classic pseudogap system [8l4l0|. There 
now exists for YBC06.5 a fairly complete, quantitative INS data set, discussed by the 
authors in terms of dynamical stripes |7(. Low-frequency data for x'KQj u )a exhibit clear-cut 
uj/T scaling, yielding the horizontal solid line in Fig. 2(a) for r e j. The INS data show a 
70/30 division between the populations of the two possible stripe domains and have a very 
nearly constant width parameter up to room temperature [7]. The YBC06.5 data differ 
from LSCO in that the (INS) values of K n j lead, through Eq.(l) and (2) with measured 
HF constants [30]], to the widely disparate dash-dot curves for r eP in Fig. 2(a). Such a 
discrepancy in the region below T c ~ 62K suggests a sharp difference in the peak widths 
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for ■x" I {q,uj) a between the INS and NMR samples. Indeed, broadening the Xi(Q,u) a peaks 
by a factor ~ 2.5 leads to unification of the r e // curves at T < 62K, as shown by the 
blue square and red circle points in Fig. 2(a). This surprising broadening effect represents 
the difference between oriented powder samples used for the 7\ measurements 2l| and 
single crystals used for the INS studies, most likely due to different degrees of chain oxygen 
ordering. Calculations of the coefficients K n i a ^ made for Fig. 2(a) used a form factor with 
discommensurations only along the a axis based on the form given by Stock et al. 7]. The 
resulting temperature-independent correlation coefficients are Kua = —0.83; Kub = —0.85; 
K 2I = 0.74; K 3Ia = 0.57; K 3Ib = 0.61. 

Values of r e ff deduced from Eqs.(l) and (2) agree very nicely with r e j data (solid line) 
in Fig. 2(a) at T c ~ 62K, again with no adjustable parameters. As with LSCO, the r e // 
curves show a sharp increase over r e i and a strong divergence from one another at T > 
62K. To find consistent values of r eP for YBC06.5, we again model x"p(<l) UJ ) using the 
squared Lorentzian form as for LSCO with the same exponential form for K\p(T). The 
result (Fig. 2(b)) is quite successful. Curves for the K n p are shown in the inset. In this case 
K\p{T) = -0.87exp[— (T — 62)/725] decays a bit more slowly and begins with a somewhat 
narrower peak (q w ~ 0.44 at T = 62K). The incommensurability ~ 0.38, however, is less than 
q w , so that the progressively broadening peak of Xp(q> u ) wn l form something of an elevated 
baseline for the incommensurate peaks. Such a background will be difficult to detect with 
unpolarized neutrons. 

The r e p curves in Fig. 1(b) and 2(b) are qualitatively similar, with r e p vanishing nearly 
linearly as T declines toward T$t r , while bending over towards room temperature. At the 
latter point, the new term contributes far more to 1/Ti than the incommensurate 'stripe' 
fluctuations. Values of r e j, which are considerably larger for LSCO than for YBC06.5, obey 
oj/T scaling and also do not display a spin gap until the materials become bulk superconduc- 
tors at T c . The experimental conclusion is therefore clear - the much-celebrated magnetic 
pseudogaps in these systems are gaps in the new term deduced from T% data, which accounts 
for more spectral weight than the incommensurate spin fluctuations at NMR frequencies. 
Moreover, given the strong evidence that the stripe (incommensurate) fluctuations which 
dominate low- and medium energy neutron measurements compete with superconductivity, 
it is the pseudogap terms that are much more likely to form a pair binding texture. 

Some time ago, in the first quantitative test of the magnetic fluctuation-dissipation the- 
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orem a joint analysis of NMR-INS data on LSCO was presented 20] . That work was only 
a partial success because of its rather simple treatment of the 7\ process. At the time, 
a two-band model was called for, but there is no longer any clear motivation for such a 



model 



311 ] . However, NMR shift analyses have been put forth recently giving evidence for 



a 'two-component' shift structure [32]. We suggest that the two-part structure for nuclear 
relaxation described in the present work could, via the Kramers-Kronig relation, form the 
basis in principle for a two-component NMR shift. In practice, there are no x"(q,uj) data 
near q = to provide a quantitative basis for a shift estimate. However, the proposed NMR 
shift structure is regarded as a natural extension of the present two-component model for 
x"(q,oj). We emphasize that the latter model does not imply two independent bands of 
charge carriers. 

In closing, we comment on the prospects for resolving the new fluctuation term x'p{q^) 
with INS methods. The squared Lorentzian model form factor used here yields estimates 
for x'piQi 1 ^) that are generally much broader and flatter than the incommensurate peaks 
reported to date. For LSCO, such a model yields an amplitude ~10% of x'li^^a for a 
scan through adjacent peaks. Such a result is compatible with experimental spectra for 
50K <T < lOOi^]. Thus, resolving Xp(<L w ) will require polarized neutrons with a far 
higher signal-to-noise ratio than what has been reported to date. Recent improvements 
in neutron scattering technique may render this feasible. However, we conclude that the 
Xp(q, term, whatever its precise shape and behavior may be, must exist, rendering x"(q, ui) 
consistent with pseudogap studies using other probes. We suggest this inference to extend 
to other cuprates as well, offering a clear and broadly consistent picture of the pseudogap 
effect throughout this family of superconductors. 

We thank B. S. Shastry and J. Zaanen for highly informative conversations. This 
manuscript has been authored by UT-Battelle, LLC, under Contract No. DE-AC05- 
00OR22725 with the U.S. Department of Energy. 
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FIG. 1: (a) Values of T e i(T) and r e // determined in various ways for LSCO are plotted as functions 
of temperature. Values of T e i(T) (solid triangles) were calculated with Eq.(4). The solid line is a 
linear regression showing approximate uj/T scaling for T > T$tr ~ 50K. Values of r e ff obtained 
from T\ data with Eq.(l) and (2) are also plotted for 63 Cu (squares) and 1 O (circles), respectively. 
The Korringa-like behavior of the r's below T$t r is highlighted by the solid red line. Data for 17 O 
are scaled to that line for the analysis in part (b). For a general comparison, values of T e ff(T) for 



YBC07 obtained by Uldry and Meier 



22] for that system are shown as a dashed line, (b) Values 



of r e p(T) obtained with Eqs.(5) and (6) using the same T\ data as above are plotted against 
temperature. For this purpose, values of the i^ n p's derived from a squared-Lorentzian model are 
used. The K n p } s with their fitted temperature dependence are shown in the inset. The solid line 
representing the behavior of r e / from part (a) is replotted here for comparison. See text for further 
discussion. 
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FIG. 2: Plots of T e jj, r e i, and r e p are presented for YBC06.5 similar to the LSCO case in Fig.l. 
(a) Values of T e j(T) for YBC06.5, calculated with Eq.(4) using INS data from Ref. 7, are shown 
as a solid line that obeys u/T scaling. Values of T e ff obtained with Eq.(l) for the 63 Cu and with 
Eq.(2) for the 17 O from Ti data [21] using values of K n j are plotted as dash-dot lines. The disparity 
between the dash-dot lines for 63 Cu and 17 is attributed to a disparity in incommensurate peak 
widths between NMR and INS samples and is corrected using adjusted peak widths (see text), 
leading to the curves showing filled squares ( 63 Cu) and circles ( 17 0). At T = 62K the latter results 

1 as agreement with r e j. (b) Values of r e p obtained with 



show good mutual correspondence as we 
Eq.(5) and (6) using T\ data from Ref.[2l| are plotted against temperature. Calculation of the 
K n p^s is described in the text. A solid line representing the behavior of T e j from INS data (see 
part (a)) is replotted here for comparison. The K n p's with their fitted temperature dependences 
are shown in the inset. 
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